IASSNS-HEP-92/46 

CfPA-TH-92/17 

NSF-ITP-92-120 

UMAC92-1 

June 1992 



(N 

< 



Instability and Subsequent Evolution of Electroweak Bubbles 



Marc Kamionkowski*^ 



^ ■ School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 

(N ■ 

O ■ and 
oo . 

O . Center for Particle Astrophysics, University of California, Berkeley, CA 94720 



and 



^ . Katherine Freese ' 



Department of Physics, University of Michigan, Ann Arbor, MI 48109 

and 

Institute for Theoretical Physics, University 
of California, Santa Barbara, CA 93109 



ABSTRACT 



* e-mail: kaniion@guinness.ias.edu. 
t SSC Fellow 

I e-mail: freese@umiphys.bitnet 

§ Presidential Young Investigator and Sloan Foundation Fellow 



Bubbles in a first-order electroweak phase transition are nucleated with radii 
Rq and expand with velocity v. If v is subsonic, a bubble becomes unstable to 
non-spherical perturbations when its radius is roughly 10^it!o- These perturba- 
tions accelerate the transition, and the effective velocity of bubble growth rapidly 
becomes supersonic. The transition should subsequently proceed spherically via 
detonation. If for some reason the onset of detonation is postponed, the surface 
area of the bubbles may be enhanced by 10^. We discuss consequences for elec- 
troweak baryogenesis. 



There has been much interest in the dynamics of a possible first-order elec- 
troweak phase transition (EWPT) recently. The motivation is clear: one of the 
fundamental problems in particle physics and cosmology is the origin of the baryon 
asymmetry of the Universe (BAU). While C and CP violation [1] as well as baryon- 
number violating instanton effects [2] have been known to exist in the Standard 
Model for quite some time, only recently has it been suggested that the rate for 
baryon-number violating interactions may become appreciable at high temper- 
atures [3]. If the EWPT is first order, the third of Sakharov's [4] criteria for 
baryogenesis, out-of-equilibrium processes, may also be found in electroweak (EW) 
physics. Thus, the BAU may be explained in the SSC era. 

In the standard picture of a first order EWPT, spherical bubbles are nucleated 
with microphysical radii Rq ~ 10~^^ cm and then expand spherically with velocity 
V to macroscopic radii Rperc ~ IO^'^ vRq before they collide. The bubble-wall 
velocity v is still uncertain, but recent estimates suggest that a wall may propagate 
subsonically {i.e., as a deflagration front) [5]. In this Letter we show that shape 
instability of the bubble wall rapidly causes the propagation to become turbulent 
and proceed more quickly, and probably instigates the onset of detonation; then 
the bubbles expand spherically to fill space at supersonic velocities. 

A deflagration front is unstable to perturbations with wavelengths in the range 
^ ^ ^ ^max where Ac ~ Rqv~'^ is set by the surface tension in the wall [6,7,8], 
and Xmax-i which is proportional to the radius it! of the bubble, is set by the un- 
derlying expansion of the bubble [7] . After the bubbles are nucleated, they expand 
spherically until Xmax reaches Ac; then Rinst ~ 100i'~^i?o and hydrodynamic in- 
stabilities set in. The subsequent bubble shape will be roughly spherical; however, 
instead of a smooth surface, the wall is highly wrinkled with distortions that en- 
hance the surface area of the wall and thereby accelerate the transition. Although 
the details are far from understood (in any fluid dynamic system), the onset of 
turbulence and corresponding acceleration of the transition should result in a det- 
onation front shortly after the bubbles are nucleated, when roughly a fraction lO"^-*^ 
of the Universe has been converted to the low-temperature phase. On the other 
hand, if for some reason the transition continues as a deflagration, then the surface 
area of the walls is enhanced by five orders of magnitude by the time the bubbles 
percolate. In either case, the dynamics of the transition assumed in models where 
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baryon number is produced in a first-order EWPT could be significantly altered. 

We limit ourselves to the case where the transition occurs at temperatures 
near the critical temperature Tc, the latent heat is small compared to the thermal 
energy density, and f ^ 1; under these assumptions the calculations simplify 
considerably. Such conditions are possible in the Standard Model [5]. We feel 
that a more general analysis should result in similar conclusions. We neglect the 
expansion of the Universe, since the timescale for the EWPT is much smaller than 
the expansion timescale [5,9]. 

First we review some results from the theory of combustion of relativistic fluids 
[10] in the case where fluid velocities are nonrelativistic. Consider a planar interface 
in the y-z plane that propagates in the —x direction. Then in the rest frame of 
the wall, matter in the symmetric phase enters the interface with a velocity vi, 
and matter in the broken-symmetry phase leaves the interface with velocity V2- 
Conservation of energy and momentum across the interface leads to the conditions 

wivi = W2V2, and pi=P2, (1) 

where w — e + p is the enthalpy density, e is the energy density, and p is the 
pressure. Throughout, the subscript "1" refers to the symmetric phase and the 
subscript "2" refers to the broken-symmetry phase. 

The ei and pi may be obtained from finite-temperature field theory. Near Tc 
the effective potential may be written [5,9,11], 

y (0, T) = (i/2)7(r2 - r2)^2 _ (1/3)^7^^3 ^ (1/4) a04^ (2) 

where is the Higgs field, Tq = {'y/X)(f)Q, (po = 250 GeV is the Higgs vacuum 
expectation value, and a, A, and 7 are parameters that depend on the W, Z, and 
top-quark masses, and on the Higgs structure of the theory. The critical temper- 
ature Tc is defined as the temperature at which there exists a second minimum of 
V degenerate with the minimum at the origin; for the effective potential above, 
T!-Ti/{l-lf,) [9]. 

Although the difference in free energies B{T) between the two phases is in 
general a complicated function of T, if the transition occurs near Tc then B(T) ~ 
(L/4)[l - (T^/T^)], where L = -Tc{dB/dT)\T^ = {Aa'^^/9X^)T^T^ is the latent 
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heat of the transition [9] . This leads to the rather simple equation of state (which 
mimics the QCD bag model), 



pi(r) = K(r)-L]/4, and ei(r) = [3«;i(r) + L]/4, (3) 

P2{T)^W2{T)/A, and e2(T) = 3^/;2(T)/4, (4) 

where wi{T) = aiT^, and W2{T) = a2T^-, oi ~ 02 ~ 100 and ai — 02 = L{4T^)~^ 
(see Ref. 9). 

We now study the hydrodynamic stability of a spherical bubble to small non- 
spherical perturbations. For distortions with A <^ i?, it is valid to treat the wall as 
a planar interface. The stability of a planar front for combustion of a relativistic 
gas was recently discussed by Link [8]. Consider a small perturbation to the planar 
discontinuity of the form Xf = dexp{iky + ut). This disortion in the wall surface 
will be accompanied by perturbations to the velocity and pressure. If v and p are 
unperturbed quantities, then the perturbations p' and if must satisfy [8] 



p' + wCgV • = 0, and 



if + Ivp = 0, (5) 
w 



where Cg = l/^/S is the speed of sound. There are four boundary conditions on 
the discontinuity that the perturbed quantities must satisfy. The first, 

follows from Eq. (1) and includes the effects of surface tension and finite mass 
density of the wall; these two effects favor a flat surface. The boundary condition 
on {wVxY (to lowest order in vi and V2,), from Eq. (1), is 

Wl - dxf/dt) = W2 {v'x2 - dxf/dt) . (7) 

Requiring that the tangential velocities on both sides be equal leads to 

v'y-^ + vi{dxf/dy) = Vy2 + V2{dxf/dy). (8) 

We make the ansatz that the enthalpy flux across the interface is proportional to 
the net blackbody energy flux across the interface. Then the perturbed enthalpy 
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flux is [8] 

M'=(3/4)a(pl-py, (9) 

where a is a fudge factor, < a < 1, and a = 1 in case the enthalpy flux is the 
blackbody energy flux. (Our flnal results will not depend on a.) Equating Eqs. (7) 
and (9) gives us our fourth boundary condition. 

If nontrivial solutions that satisfy Eq. (5) and the boundary conditions can 
be found for some a; > 0, then there are growing modes, and the wall is unstable 
to small perturbations. To satisfy the equations of motion and the boundary 
conditions, cu must satisfy [8] 



,2/ N cr^^ 

k [vi - V2) + 



WlVl 



ViV2 = 0. (10) 



If V2 > vi {i.e., if the phase transition proceeds via deflagration), then Eq. (10) has 
a positive root for wave numbers k < kc — {v2 — vi)wivi/a, and there are growing 
modes. For larger wave numbers the system is stabilized by surface tension. 

Prom Eqs. (1), (3), and (4), we flnd V2 — vi ~ {L/wi)vi; thus the wall is 
unstable to small perturbations with 

X>Xe = k-^c^a/{Lvl). (11) 

In the limit of small supercooling the surface tension is given by cr = (2"^/^a;'^/3'^A^/^)Tg 
[9] . Determining vi is much more difficult and requires an investigation of the mi- 
croscopic interactions of the particles in the thermal bath with the advancing wall. 
Recent estimates suggest that for the minimal standard model, the wall velocity 
may be in the range 0.01 ^ vi < 0.3 [5]. Note that Ac ~ RqVi^ ~ 10~^^ cm (for 
vi ~ 0.1). 

As long as A <^ i?, the analysis assuming a planar interface should be valid; 
however, for A ~ i?, one should take the expansion of the bubble into account. As 
the bubble expands, the wavelength of the perturbation increases. If the amplitude 
of a perturbation grows more slowly than the wavelength, then the distortion is 
smoothed out in time. For our case of a weak transition, 5 = {vi—V2)/vi ~ L/w 
0.01 <^ 1, the growth rate for an instability with a large A is found from Eq. (10) 
to be a; ~ 5v\k/2, while the growth rate for the bubble is roughly vi/R. For a 
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perturbation to be to be unstable it must have k > 2/{R6), while perturbations 
with A > Xmax — RS/2 will be stabilized. Although this derivation is heuristic, 
for 5 <^ 1 it reproduces the results of the exact analysis for the case of spherical 
combustion of a nonrelativistic gas, and it should be a good approximation in the 
case of a slowly-moving relativistic gas as considered here. So, perturbations with 
A < Ac will be stabilized by surface tension, and those with A > Xmax will be 
stabilized by the growth of the bubble. 

In the standard picture of a first-order EWPT, bubbles are nucleated with radii 
Rq cr/L and then grow spherically with velocities v until the bubbles percolate. 
The radius of the bubble at this time is (again, assuming typical parameters [9]) 
Rperc ~ ^0^^ v{Mpi/Tc)T^^ ~ IO^^vRq. On the other hand, when Xmax reaches 
Ac, Rinst — crw / {L^v^) ~ 100f^^i?o; at this point the bubble becomes unstable to 
non-spherical perturbations. Since Rperc is niany orders of magnitude larger than 
Rinst the perturbations have plenty of time to mature and the standard picture of 
bubble evolution may be drastically altered. 

If the bubble volume is V ^ the bubble will look somewhat spherical with a 
nominal radius R given by V = (4/3)7ri?^; however, instead of a smooth surface, 
the wall is highly wrinkled on scales Ac ^ A < Xmax, and the surface area of such a 
bubble is actually much larger than AuR^. Perturbations to the fiuid fiow accom- 
pany those in the bubble surface, so that the normal fiow velocity of fiuid across 
the interface is vi at every point on the surface [c/. Eqs. (7) and (9)]; therefore, the 
rate of the transition is enhanced. A similar situation arises in supernova theory, 
where the burning of a carbon-oxygen white dwarf proceeds via defiagration and 
the rate at which burning occurs is proportional to the surface area of the wrinkled 
flame [12]. The surface area is enhanced roughly by a factor 



where the fractal dimension D [13] is some number between 2 and 3 but most 



Although the exact fractal dimension is uncertain, the qualitative form of Eq. (13) 




(12) 



hkely near 2.6 [12]. The effective velocity v^ff = (dR/dt) at which a sphere of 
comparable volume as the bubble would expand becomes 

Veff ^ {Xmax/Xc)^~'^ V2. (13) 
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is correct; Vf,ff might differ from our estimate by an order of magnitude or so, but 
tliis lias little effect on our conclusions. There is a possibility that once the pertur- 
bation goes nonlinear {i.e., its amplitude becomes comparable to its wavelength) 
that it becomes stabilized and the resulting flow is not turbulent [7]. However, in 
this case, the wall would still be wrinkly on length scales from to Xmax- The 
resulting surface-area enhancement and v^ff would still be comparable to those 
given in Eqs. (12) and (13). Shortly after instabilities set in, the transition acceler- 
ates and when R ^ ~ Rqv^^6^^ <^ Rperci Vf, ff becomes supersonic. 
At this point only a fraction (R/ Rperc)^ ~ 10~^^ (for fi ~ 0.1) of the Universe has 
been converted to the new phase. Therefore, if baryogenesis occurs at the EWPT, 
it takes place after Wgjj becomes supersonic. 

The most likely scenario is that when Vgff increases past Cg, a detonation wave 
sets in. Simply stated, the reason is that the deflagration front is preceded by a fluid 
flow, and it is hard to see how the appropriate fluid flow can be maintained in front 
of a deflagration wave itself moving supersonically. In the frame of the deflagration 
front, the flow velocity of fluid into the interface is smaller than the flow velocity of 
fluid out of the interface {i.e.,V2 > vi). In the "rest" frame of the Universe, the fluid 
is at rest far away from the transition; furthermore, by symmetry arguments, the 
fluid inside the bubble must be at rest. There is a piston effect as the wall pushes 
the fluid outside the bubble with a speed V2 — vi — 5vi. A precompression shock 
precedes the deflagration front and accelerates the fluid, which is initially at rest, 
radially outward to a velocity Svi [10]. If the wall is distorted and the transition is 
accelerated, the wall pushes the fluid outside the bubble with a velocity near Svgff. 
Since Svgff <^ 1, only a weak shock is needed. Weak shocks travel at velocities 
only slightly larger than Cg [14]; once v^ff > Cg, the deflagration front will merge 
with the shock to form a detonation wave. Although the exact mechanism for 
onset of detonation from deflagration is still under investigation in fluid systems 
and is not entirely understood even in the nonrelativistic case [7,15], the onset 
of detonation from a shock preceding an accelerating turbulent deflagration front 
is observed in laboratory experiments [7], and appears in the theory of Type la 
supernovae [15]. In order to satisfy the hydrodynamical equations of motion with 
the boundary conditions that the fluid far from the bubble as well as at the center 
of the bubble be at rest, the Chapman-Jouget condition must be satisfled, and the 
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bubble expands at a velocity v c:=^ ^1/3 + \/25/9 (for 5 <^ 1) slightly larger than 
Cs [14]. Since perturbations cannot propagate faster than Cg, perturbations should 
be smoothed out and the bubble expands spherically to fill all space. 

A detonation wave heats the gas as it passes, so one might worry that if the 
gas is heated to a temperature above Tc that the phase transition cannot continue. 
A detonation will certainly propagate if the supercooling of the Universe is greater 
than the subsequent heating, as is found in some models (though not all). Further 
work should investigate the details of the phase transition at a detonation front in 
the case that the gas is heated above Tg. 

On the other hand, if for some unforeseen reason, a terminal v^ff smaller 
than Cs is reached, then the subsequent evolution could continue as a deflagration 
with a distorted surface. This distortion can be quite dramatic: By the time 
the bubbles percolate, the surface area of the bubbles is enhanced by roughly 
{RpercS / Xc)^~'^ ~ 10^ [c/., Eq. (12)]. If this is the case it might play a role in EW 
baryogenesis. However, the baryon number in recently proposed existing models 
where baryogenesis occurs at the phase boundary [16,17] should be unaltered. 
Although the transition in this case would be accelerated, the resultant baryon 
number is generally proportional to the amount of fluid that passes through the 
wall and this remains unaltered by turbulence. In some models such as that in 
Ref. 17, the rate of baryogenesis depends on transport of particles near the wall 
and the flnal baryon number depends on the wall velocity; in such models, the 
resulting baryon asymmetry depends not on v^f f but only on the local velocity of 
fluid through the wall, which remains unchanged (in the nonrelativistic limit) by 
turbulence. We can only speculate that relativistic corrections could actually alter 
the flow velocity across the wall. Another possible effect of the wall convolution 
is that in models where transport near the wall is crucial, particles could multiple 
scatter off one wall into another wall; however, this would require that Ac be smaller 
than the particle mean-free path. 

Throughout we assumed that latent heat is transported from the surface hy- 
drodynamically. If, on the other hand, radiative transport is important and bubble 
growth is limited by diffusion of latent heat from the wall, then the wall may be- 
come unstable on length scales larger than mean-free path of radiation as shown 
by Preese and Adams [18] for the case of a flrst-order QCD phase transition. If 
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so, instabilities may set in even earlier than we found (as soon as the bubbles nu- 
cleates), and the surface- area enhancement could possibly be even larger than we 
estimated. The resulting bubble shape in this case may deviate drastically from 
spherical; the bubble looks like seaweed. Multiple scatter would also become more 
important in this case. 

To summarize, the propagation of a deflagration front in a weakly first order 
EWPT becomes turbulent, the transition is accelerated, and the effective propa- 
gation velocity of the walls rapidly becomes supersonic. Under these conditions 
the deflagration front could turn into a detonation shortly after the bubble is 
nucleated, and the macroscopic growth of the bubbles should occur via a deto- 
nation wave traveling near Cg. Our results suggest that, due to hydrodynamic 
effects, macroscopic bubble propagation may be significantly different from what 
one would expect from detailed studies of the microscopic kinetics [5]. This should 
come as no surprise; it has long been known that the propagation velocity of a 
spherical detonation wave is determined by hydrodynamics (the Chapman- Jouget 
condition [6,10,14]) and not by the microscopic kinetics of the reaction. 

Strictly speaking, our analysis is valid only for nonrelativistic propagation ve- 
locities and for transitions with small latent heat, but a more general analysis under 
less restrictive assumptions should result in similar conclusions. For example, as 
long as the deflagration velocity v is subsonic, small perturbations could propagate 
ahead of the detonation front, and the hydrodynamic instability should exist; in 
addition, as v is increased (while still subsonic), is decreased [c/., Eq. (11)] so 
the instability should set in sooner. If the latent heat is increased, Ac decreases and 
Amaa; becomes larger, so turbulence should set in sooner. Also, for larger latent 
heats, the detonation front propagates at a larger velocity [14]. The transition from 
deflagration to detonation may also be important for the dynamics of the QCD 
phase transition if it is first order. 
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